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On the example of the three-dimensional Ising model, we show that nonperturbative renormaliza- 
tion group equations allow one to obtain very accurate critical exponents. Implementing the order 
d 4 of the derivative expansion leads to v = 0.632 and to an anomalous dimension r\ — 0.033 which 
is significantly improved compared with lower orders calculations. 
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Many problems in high-energy as well as in statistical 
physics call for nonperturbative methods. On the one 
hand, several physical systems are described by field the- 
ories in their strong coupling regime so that the usual 
perturbative techniques become troublesome. They fail 
either because only the first orders of perturbation are 
computed and do not suffice, or because, even when high 
orders are known, standard resummation techniques do 
not provide converged results. On the other hand, some 
phenomena such as confinement in QCD or phase tran- 
sitions induced by topological defects are genuinely non- 
perturbative. 

Apart from some methods restricted to specific di- 
mensions or situations, very few nonperturbative tech- 
niques are available. During the last years, the Wilson 
approach 1 to the renormalization group (RG) has been 
turned into an efficient tool 23 4 . This nonperturbative 
RG can be implemented in very general situations and, 
in particular, in any dimension, so that it has allowed 
one to study several issues difficult to tackle within a per- 
turbative framework among which the three-dimensional 
Gross-Neveu model 5 -, frustrated magnets 6 -, the randomly 
dilute Ising modeS, and the Abelian Higgs model^. 

This method relies on a nonperturbative renormaliza- 
tion of the effective action T, i.e. the Gibbs free energy. 
It consists in building an effective action Tk at the run- 
ning scale k by integrating out only fluctuations greater 
than k. At the scale k = A, A -1 denoting the spacing 
of the underlying lattice, Tk coincides with the Hamil- 
tonian H since no fluctuation has yet been taken into 
account while, at k = 0, it coincides with the standard 
effective action T since all fluctuations have been inte- 
grated out. Thus, Tk continuously interpolates between 
the microscopic Hamiltonian H and the free energy T. 
The running effective action Tk follows an exact equa- 
tion which controls its evolution with the running scale 
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where t = ln(fc/A) and 1^ '[<p\ is the second functional 
derivative of Tk with respect to the field <f>(q). In Eq. QJ, 



Rk (q) is an infrared cutoff function which suppresses the 
propagation of the low-energy modes without affecting 
the high-energy ones. 

Although exact, Eq. JJJ is a functional partial integro- 
differential equation which cannot be solved exactly. To 
handle it, one has to truncate Tk- A natural and widely 
used truncation is the derivative expansion, which con- 
sists in expanding Tk in powers of d(j>, keeping only the 
lowest order terms. Physically, this truncation rests on 
the assumption that the long-distance physics of a given 
model is well described by the lowest derivative terms, 
the higher ones corresponding to less relevant operators. 
Up to now, only truncations up to order d 2 have been 
considered since, in many cases, they turn out to be suf- 
ficient to get a satisfying qualitative and even sometimes 
quantitative description of both universal and nonuniver- 
sal behaviors^. 

Nevertheless, several important issues concerning the 
reliability of the method remain open. The first one con- 
cerns the convergence of the derivative expansion. This 
point is particularly delicate since, within this kind of 
truncation, there is no expansion parameter in terms 
of which the series obtained can be analyzed and con- 
trolled. It has been moreover suggested that the expan- 
sion could be only asymptotioifi. Actually, this question 
of convergence has only been addre ssed within the per- 
turbative context. In Refs. |llll2fl . it has indeed been 
shown that the two-loop perturbative result for the O(N) 
model can be recovered from a summation of the deriva- 
tive expansion. However, in its full generality, this prob- 
lem still appears as a major challenge. The second is- 
sue, the truly interesting one from a practical point of 
view, concerns the accuracy of the results provided by 
low orders truncations. This has been thoroughly stud- 
ied only at order d° within the O(N) model and at order 
d 2 within the Ising model through the optimization of 
the cutoff function Rjjl&ll (see also related studies with 
the Polchinski equatio n 15 ' 16 and within the proper time 
RG formalismii) . Let us emphasize that, even for these 
models, the anomalous dimension r\ remains poorly de- 
termined. This likely originates in the crudeness of the 
order d 2 truncation that fails to capture the essential 
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momentum dependence of the two-point correlation func- 
tion. In this respect, an important remark is that in the 
critical theory, and at k = 0, this function is nonanalytic 
— rj^2 Q (g) ~ q 2 ~ v — , so it appears nontrivial to retrieve 
77 from a derivative expansion. However, for k 7^ 0, the 
infrared fluctuations (q <C k) are suppressed and T K k '(q) 
should become regular with the standard q 2 behavior. 
This means that the nonanalyticity builds up smoothly 
as k vanishes. Berges et al£ have proposed that this func- 
tion behaves approximately as q 2 (q 2 +ck 2 )~ v ^ 2 , where c is 

(2) 

a constant. Roughly speaking, for T y k ' , the derivative ex- 
pansion consists in expanding this function around q = 
and in computing 77 from its behavior in fc 2 , instead of 
q 2 . It is not trivial that the resulting series for 77 con- 
verge since it amounts to correct the normal q 2 behavior 
with higher powers of q 2 . The aim of this paper is to 
investigate this question by including order d terms in 
the derivative expansion of Tk for the three-dimensional 
Ising model. 

The effective average action Tk of the Ising model trun- 
cated at order <9 4 is written as: 



d d x { U k (p) + 



z k ( P ) (y<py + w%(p) (A<t>y 
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w b k { P ) (v^r(^)+w^p) (V0) 



where p = (f> 2 /2 is the Z2 invariant. Compared with 
its expansion at order d 2 , T k involves three new terms 
denoted W k (p) s = a,b,c, linearly independent with re- 
spect to the integration by parts. The evolution equation 
for the potential U k is derived by evaluating Eq. (JIJ for 
a uniform field configuration. By contrast, the definition 
and, thus, the evolution of the functions Z k and Wj!'s are 
linked to a specific momentum dependence of the func- 
tional derivatives of T k , in the limit of vanishing external 
momenta: 

<5 2 r 

Zk{p) = ^ mfttov (3) 
wa ^ = ^jmk^y 3 (4) 

wbk{p) = -^A^wSww' (5) 



(6) 

As usual, to find a fixed point, we use the associated 
dimensionless renormalized quantities p,Uk, z k and wf.. 
The flow equations of these functions, derived from Eq. 
JTJ, are far too long to be displayed. 

As in Ref . |14| , we have implemented a further approx- 
imation which consists in expanding each running func- 
tion Uk, z k and the w^s in powers of p. The motivation 



which underlies this is twofold. First, in systems hav- 
ing a symmetry group smaller than 0(A), the number 
of functions analogous to u k , Zk and wf. grows as well as 
the number of arguments, analogous to p, on which they 
depend. In this case, dealing with the full field depen- 
dence at each order of the derivative expansion can be 
very demanding and the field expansion becomes almost 
unavoidable. Second, this expansion provides valuable 
indications about the orders in field necessary to correctly 
describe the critical behavior. Here, we expand the u k , 
Zk and functions around the configuration p — p that 
minimizes Uk since it leads to a better convergence than 
the expansion around p = 



(a. 



Pi 



(7) 



where £ stands for u, z, w a , w b , w c . The RG equation 
JU then leads to a set of ordinary coupled differential 
equations for the coupling constants {Cj.fc}- The nonper- 
turbative features of their flows with the running scale 
k are entirely encoded in a finite set of integrals, called 
threshold functions. There are six - three of them being 
specifically linked to the inclusion of the d A order terms 
- which are written: 



F*= / dyyz^dt (f(y) 



1 



(p{y) + m 2 ) n 



(8) 



where dt means that the derivative only acts on the cutoff 
function R k {q) = Z Q<k q 2 r(y) with y = q 2 /k 2 ; p(y) = 
2/(1 + w% k y + r(y)), m 2 = 2u 2 ,k pa, and f(y) can be 
either y(d y p) 1 with i = 0, ...,4 or yd 2 p. The occurence 
of d 2 p imposes the cutoff function R k to be at least of 
class C 3 , which dismisses for instance the theta cutoff 
introduced ini£. Here, we choose the exponential cutoff 
defined b}^: 



Rk{q) = a- 



Zp, k q 2 



(9) 



which fulfills this condition and constitutes an efficient 
regulator. We remind that any truncation of Tk intro- 
duces a spurious dependence of the results on R k . Here, 
we study this influence by varying the cut-off through the 
amplitude parameter a. For each truncation, the optimal 
a is determined through a principle of minimum sensitiv- 
ity (PMS) which indeed corresponds to an optimization 
of the accuracy of the critical exponents^. 

At each order of the derivative expansion, up to or- 
der d A , and for higher and higher order field truncations, 
we compute the fixed point and the associated critical 
exponents v and r\, as functions of a. Then, for each 
truncation, we determine the optimized exponents from 
the PMS, which are referred to, in the following, as PMS 
exponents. We first expand in fields the potential Uk, and 
then Zk, which respectively constitute the orders d° and 
d 2 of the derivative expansion. The corresponding PMS 
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FIG. 1: i/pnis and ?7pms as functions of the truncation. The 
three zones I, II and III correspond to the expansions of Uh, 
Zk and the w%"s respectively. In zone III, the two values at 
p w s = 3 reflect the different choices of PMS solutions (see 
below) . 



exponents are displayed as functions of p u>z — which de- 
note the orders of the truncation in p of Uk and Zk — in 
the first two zones of Fig. 2] At this stage, it is worth em- 
phasizing that strong oscillations occur at the first orders 
in the field expansion for both orders d° and d 2 . It fol- 
lows that the PMS exponents become almost steady only 
from p U;Z — 4. As discussed in Ref. 14], the truncation 
Pu = 8 and p z — 6 allows one to obtain a very accurate 
approximation of the order d 2 results. Indeed, the cor- 
responding exponents i/ PMS = 0.6291 and t^pms — 0.0440 
differ by less than 1% compared with their "asymptotic" 
values obtained for large p UjZ (see Table©. Note also 
that, already at this order d 2 , v FMS agrees well with the 
best known values whereas, as mentioned above, this is 
not the case for 77 PM s- 

Let us come to the role of the order <9 4 terms. We 
choose to simultaneously expand in fields the three func- 
tions Wfc, s = a,b,c, up to p w s = 5, while fixing p u = 8 
andp z = 6. Actually, the highest truncation corresponds 
top„,c = 5 and p w a.b = 4 for the following reason. Figure 
|5]displays, for each wf, considered independently, the evo- 
lutions of the PMS exponents with the order of the field 
truncation. It shows that the exponents associated with 
or w\ have almost converged, up to a few tenths of 
percent, as soon as = 3. On the contrary, ?7pms and 
i/ PMS related to still oscillate at this order. We have 
checked that, within the simultaneous expansion of the 
three wt, indeed dictates the variations of the critical 
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FIG. 2: z/pms and ?7pms as functions of the order of the field 
truncation for each function wl separately. For w c , the two 
values at p w c = 3 reflect the different choices of PMS solutions 
(see below). 



exponents, mi°j. and Wj k exerting a minor influence for 
3 

This, together with the fact that we encounter here 
the limits of our computational capacities, justifies our 
choice (pw a ,Pw b iPw a ) — (4,4,5) for the last truncation. 

We can now concentrate on the behavior of the expo- 
nents at the order d A . At the low orders field trunca- 
tions, corresponding to p w $ = 0, 1, and 2, each exponent 
exhibits a single PMS solution, ^ PMS and ?7pms, which 
are thus unambiguously defined. As displayed in Fig. [3] 
several PMS solutions appear for the next two trunca- 
tions, corresponding to p w s — 3, 4. This renders the op- 
timization procedure in these cases (see discussion below) 
unclear. Concerning the largest truncation, 7y PMS is un- 
ambiguously determined from the unique PMS solution. 
For v, several PMS solutions exist. However, provided 
the field expansion has almost converged at this order, 
a unique PMS solution can also be selected for v. The 
argument underlying this choice originates from the fact 
that when no truncation in derivatives is performed, the 
results are independent of the cutoff. Therefore, the best 
cutoff is the one achieving the weakest sensitivity of the 
results with respect to the order of the derivative expan- 
sion, i.e. leading to the fastest convergence (see Ref. 
for a detailed discussion). In practice, this consists in 
minimizing the difference between the values of v deter- 
mined at order d n and at order d n+1 . In our case, this 
selects the PMS solution located at a ~ 0.6 (see Fig. [3J. 

Let us now discuss the convergence of the field expan- 
sion. To this end, we first examine the two truncations 
p w s = 3, 4 for which multiple PMS solutions exist for 
both exponents. There is no argument to clearly settle 
between the PMS solutions. We present two sensible way 
to favor reasonable PMS solutions. First, one can choose 
to minimize, for these orders, the oscillations induced 
by the field expansion. This, in turn, corresponds to im- 
proving the rapidity of convergence of the field expansion. 
This choice corresponds to the full line in the third zone 
of Fig. n Alternatively, one can decide to follow a given 
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FIG. 3: Variations of the critical exponents v and r\ with the 
cutoff parameter a, Eq. for the highest orders truncations 
labelled by (p-w a ,p w b,pw c )' 



ously steady, this suggests that the asymptotic regime 
is just entered. This is consistent with the fact that, 
at orders d° and <9 2 , the oscillations die down for the 
same order of truncation: p u _ z ~ 4. This legitimates our 
former assumption of field convergence. We therefore ap- 
proximate the order d 4 results by the p w s = 5 estimates, 
see Table I. To summarize, we have computed the criti- 
cal exponents of the three-dimensional Ising model up to 
the d 4 order in the derivative expansion. The successive 
contributions significantly decrease with the order, which 
supports good convergence properties of this expansion, 
and in particular a correct behavior of its implementation 
around q = 0. We emphasize that the exponent v is al- 
most unaltered at order d 4 compared with its value at the 
order d 2 , whereas rj undergoes a substantial correction 
which drives it within a few percents of the best known 
values. This confirms the statement that the inclusion of 
the d 4 order terms allow one to improve the anomalous 
dimension. Note that although fully converged results 
would require to handle the full field dependence of Uk, 
Zk and the w s k , fM> 1 this study shows that the truncation 
in fields constitutes a reliable way to compute critical 
exponents. Finally, the present work brings out convinc- 
ing evidence of the ability of the effective average action 
method to provide very accurate estimates of physical 
quantities. 



PMS solution (characterized by its concavity and its lo- 
cation a), order by order in the field expansion, starting 
from the low orders p w > = 0, 1,2, where the PMS solu- 
tions are unique. This corresponds to the dashed line in 
Fig. ^ Note that both criteria lead to the same PMS so- 
lution for p w s = 4. Finally, the important features of the 
exponents evolution remain essentially unchanged what- 
ever choice is adopted: low orders generate strong oscil- 
lations that tend to vanish after a few orders. Indeed, the 
results for p w > = 4: i/ PMS = 0.6234 and 77 PM s = 0.0289, 
are very close to those for p w s = 5: ^ PM g = 0.6321 and 
?7pms — 0.0330. Although the exponents are not rigor- 



method 


V 


V 


LPA (a) 


0.6506 





Q2 (a) 


0.6281 


0.0443 


Qi (b) 


0.632 


0.033 


7-loop (c) 


0.6304(13) 


0.0335(25) 


MC (d) 


0.6297(5) 


0.0362(8) 



TABLE I: Critical exponents of the three-dimensional Ising 
model: a) effective average action method (field expansion) 14 ; 
b) present work; c) 7-loop calculations—; d) Monte-Carlo 
simulations—. 
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